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I. INTRODUCTION 

The standard regularization and renormalization 
procedures are very well documented for any pertur- 
bative calculation of physical observables. They rely 
however on identification and cancellation of diver- 
gences by appropriate counterterms. While these pro- 
cedures can be done in principle to any order in per- 
turbation theory, it becomes particularly difficult, and 
cumbersome, to implement them beyond leading or- 
ders. 

We shall advocate in the following study the use of a 
general regularization/renormalization scheme based 
on the intrinsic properties of quantum field opera- 
tors. Indeed, the fundamental objects of relativistic 
field theory are covariant field operators recognized 
long ago as operator valued distributions (OPVD). 
These distributions are defined on test functions with 
well-defined properties [1]. The formulation of the S- 
matrix, well documented in [21 [3j, was developed in 
the context of OPVD in the 1970s by Epstein and 
Glaser |H \S\ and more recently by Scharf [Qj . More- 
over, when causality conditions are imposed on test 
functions, S-matrix amplitudes, defined as time or- 
dered products of OPVD, are split into causal ad- 
vanced and retarded pieces, in such a way that any 
singular behavior at equal space-time points is avoided 
[7]. More recently, the connection of this approach 



with the well known BPHZ renormalization scheme 
has been shown in 0. 

We show in this study how this new regular- 
ization/renormalization scheme can be applied in 
practice to physical systems in a nonperturbative 
framework. This scheme is quite general, and applies 
to any formulation of quantum field theory. As 
we shall see in the following, it is very well suited 
when relativistic bound state systems are formulated 
within light-front dynamics (LFD) [9| and its covari- 
ant formulation (CLFD) [lOj [11]. In order to have 
a coherent presentation, we shall briefly recall the 
general properties of this approach, in the spirit of 
performing practical calculations. 

The formulation of quantum fields as OPVD 
has recently been revisited by Grange and Werner 
(GW)[TH. Requiring that the resulting field theory 
should 

i) be independent of the form of the test functions 
used for the construction of the physical fields, and 

a) preserve the basic Poincare and Lorentz invari- 
ances, 

GW are led to introduce test functions in the form 
of partitions of unity (PU) [131 [H] (see Appendix |A] 
for the definition of a PU). With this choice of test 
functions, the physical fields are mathematically well 
defined and show a familiar canonical form. 
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By the very nature of the PU construction, an ar- 
bitrary scale, characterising the building blocks of the 
decomposition of unity, comes naturally into the pic- 
ture. Such a scale is the corner stone permitting the 
use of extended distributions in rcnormalization group 
studies. 

Moreover, the test functions are chosen as super 
regular test functions (SRTFs), i.e. they are functions 
of finite extension - or finite support - vanishing with 
all their derivatives at their boundaries, either in the 
ultraviolet (UV) or in the infrared (IR) domains. The 
ultimate goal of the use of a PU-SRTF is to extend 
its support to infinity (UV case) or to zero (IR case) 
in such a way that any physical amplitude remains 
finite when the PU is equal to 1 everywhere in the 
integration domain. 

For a renormalizable field theory - e.g. with only a 
finite number of primitive divergences to be handled - 
the SRTF specification is not compulsory in principle 
since the test function, and only a finite number of 
its derivatives, should vanish at the singularities [4 . 
However, the choice of SRTFs lead to a very conve- 
nient and systematic construction of the extension of 
singular distributions in the spirit of the analysis of 
Epstein and Glaser. From the logical point of view, 
the introduction of SRTFs which are also PUs has not 
the character of an arbitrary restriction in the choice 
of test functions but is a mathematical necessity, if 
one wants to satisfy the two basic conditions (i) and 
(a) mentioned above. 

The key feature of any SRTF is its equality with 
its Taylor remainder of any order. This remain- 
der, written in terms of Lagrange's integral formu- 
las, permits in turn convenient mathematical oper- 
ations leading to the extension of singular distribu- 
tions. We shall therefore call this general regular- 
ization/renormalization scheme the Taylor-Lagrange 
renormalization scheme (TLRS). In nonperturbative 
calculations, singular contributions cannot be singled 
out a priori and their treatment is a major issue. The 
OPVD approach offers the immediate advantage of 
working with finite and mathematically well defined 
amplitudes. 

Using SRTFs, we can handle, in principle, diver- 
gences of any degree. In this sense, the TLRS satisfies 
Weinberg's statement that "the so-called nonrenor- 
malizable theories are actually just as renormalizable 
as renormalizable theories" [TS]. In both cases, all 
relevant amplitudes will be finite from the very begin- 
ning, but scale dependent. The scale independence of 
physical observables will however only be achieved if 
all relevant terms in the effective Lagrangian is con- 



sidered. The number of such terms is finite in renor- 
malizable theories, but infinite in nonrenormalisable 
ones. In a given kinematical domain, one may how- 
ever expect that only a limited number of such terms 
is dominant, leading to the predictive power of effec- 
tive field theories. In that case, Weinberg's statement 
is fully satisfied. 

Following the general properties of LFD, the state 
vector of any bound state system can be decomposed 
in Fock sectors. For practical calculations however, 
this Fock decomposition is truncated to a given order 
N, where N is the maximal number of particles in the 
Fock sectors. Each Fock sector is then described by a 
many-body wave function, called the Fock component. 
As we shall see in Sec. IIIIl the construction of these 



many-body Fock components can be extended, in a 
very natural way, in order to include the test functions 
necessary to the regularization/renormalization of the 
physical amplitudes, as mentioned above. This leads 
to a very transparent and general formulation of LFD, 
with all amplitudes being finite but dependent on an 
arbitrary scale. 

Because of this intrinsic scale, we can do a renor- 
malization group (RG) analysis of our result, provided 
our renormalization scheme is well adapted to a Fock 
state truncation, and the Fock state decomposition it- 
self is convergent. As compared to the Pauli-Villars 
(PV) regularization scheme widely used in LFD cal- 
culations, the regularization procedure inherent to the 
TLRS has several advantages. First of all, we consider 
in the construction of the state vector only physical 
states, and we do not have to enlarge the Fock space 
to deal with (negatively normed) PV fields, with one, 
two or more PV fields depending on the type of sin- 
gularities. Moreover, all amplitudes are finite from 
the very beginning, and numerical calculations should 
thus be easier to control. 

The plan of this article is the following. We recall 
in Sec. [lljthc main properties of the Taylor-Lagrange 
renormalization scheme. We then apply this scheme 
to light front dynamics in Sec. |III[ Perspectives are 
presented in Sec. |IV[ We detailed in Appendix [A] the 
construction of the test function. The explicit calcu- 
lation of the Fock components in CLFD is done in 
Appendix IS] 



II. TAYLOR-LAGRANGE 
RENORMALIZATION SCHEME 

We follow in this section very closely the develop- 
ments made in Ref. [12^. We start, for simplicity, from 
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the general solution of the Klein-Gordon equation. It 
is given by a distribution, (f), which defines a func- 
tional, <I>, with respect to a test function p [13] ac- 
cording to (in D dimensions) 

a>(p)=<0,p>= j d''v<t>{y)p{y) . (1) 

The physical field (p{x) is then defined in terms of the 
translation, Tx- of '^{p), which, in flat space, is given 

by 

^{x) = Tx^{p) = j d^2/</)(2/)p(x - y) . (2) 

Because of the properties of partial integrations with 
Schwartz test functions, f{x) obeys also the original 
field equation, Klein-Gordon in this case. 

Since the test function is a symmetric function of 
X — y, its Fourier transform writes 

Pi^-y)- J (0e'''■(^-^)/(go^q2). (3) 

The field (t>{x) can thus be decomposed, in momentum 
space, in terms of creation and destruction operators. 
After integration over pq, we get 

with = + . Clearly, to have a procedure inde- 
pendent of the choice of test functions, the convolution 
of f{x) itself according to ([2| should leave (p{x) un- 
changed [12]. We shall come back to this point later 
on. 

It is apparent that test functions should be attached 
to each fermion and boson fields, when deriving the 
effective Hamiltonian. Each propagator being the con- 
traction of two fields should be proportional to . In 
order to have a dimensionless argument for f, we shall 
introduce an arbitrary scale A to "measure" all mo- 
menta. A can be any of the masses of the constituents. 
In order to deal with massless theories, we shall take 
some arbitrary value. The final expression of any am- 
plitude should be independent of A. In CLFD, the 
test function is thus a function of only. One can 
already notice, at this stage, that there is no differ- 
ence between A and rjA, with 77 a positive number, 
and the test function should embody this important 
scaling information. 



A. Singular distributions, Taylor's remainder 
and Lagrange's formulas 

Let us assume that T{X) is a distribution with a 
singularity of order k at the origin of M , where k is 
defined by 

k = inf{s : lim^AT(AA) = 0} - D . 

To give an example, let us consider, for £> = 1, the 
distribution T{X) behaving like ^ with I — 1,2,3... 
when X ^ but with faster decrease than ^ at infin- 
ity. We have k — l + l — l — I. Any physical amplitude 
is written in a schematic way as 

r dX T{X) f{X) . (5) 
Jo 

If the test function f{X) does not tend to zero 
fast enough at the origin, the classical mathematical 
method to give a meaning to A is to perform a Taylor 
expansion of / around X = and to suppress from / 
as many terms in the series as is necessary to obtain 
a finite contribution from the region around X — 0. 
The resulting function is just the Taylor remainder at 
X = 0, denoted by R''{f), and defined by 

^ vn 

R\f){X) = f{X) ^5"/WU=o . (6) 

This method has been known in functional analysis 
for a century and leads to the concept of pseudofunc- 
tions [T3]. In the diverging functional A, f{X) is thus 
replaced by its Taylor remainder R^{f){X) of order 
k. It is also at the heart of the BPHZ renormalization 
scheme [IF. 

However, when f{X) is a SRTF vanishing at the 
origin with all its derivatives, / is strictly equal to 
its Taylor remainder of any order. In that case, the 
physical amplitude obeys the identity 

/>oo />oo 

A= dX T{X)f{X) - / dX T{X)R\f){X) . 
Jo Jo 

In order to define the extension, denoted by T{X), of 
the distribution T(X) over the entire space, i.e. in the 
limit where / ^ 1, we can use Lagrange's formula for 
the Taylor remainder, hence our denomination Taylor- 
Lagrange renormalization scheme. Their use was first 
advocated by Brunetti and Fredenhagen [T7]. The 
explicit procedure, which involves a transposition of 
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R''{f){X) to T{X), is detailed in the following. It 
leads to 



^ = lim / dX T{X) f{X) 



(8) 



The limit / ^ 1 can now be done safely since the in- 
tegral of T over X is finite. Note that the functional 
T is defined modulo a sum of the distribution 5{X) 
and its derivatives of order less or equal than k. This 
sum is killed by / in the functional, and should there- 
fore not be considered in This is however not the 
case for the expression of T{X) considered in isolation 
of the integral. A similar analysis can be done when 
X — > oo, i.e. in the UV domain. 

Lagrange's formula can take several forms depend- 
ing on which kinematical domain, UV or IR, we are 
interested in. We shall separate in the following the 
test function used in the UV domain, , from the 
one used in the IR domain, in order to separate 
their different scaling behavior. In the UV domain, 
we have for instance 



dtdt[f>{xt)\ 



(9) 



for the simplest case where the function alone is 
zero at infinity. One may rewrite this identity in a 
slightly different form, by shifting the derivative on t 
to a derivative on X. We get 



df 

■{X)^~X -dx[f>{Xt)\ . (10) 



Note that the Lagrange formula can be written, in the 
case of a function of X^ : 

f>{X')^-X j^^jdx[f>{Xh')] . (11) 

It can be generalized if the function is a SRTF, 
leading to 

f(l-«+M^V>(^t)], (12) 

for any integer fc > 0. 

Similarly, one can write in the IR domain: 

X dt 



= - J^ -(1 - t)'^d'/' [X''f<{Xt)] . (13) 

Lagrange's formula can be written in different ways 
since, if f{X) is a SRTF, then X''f{X) is also a SRTF 



for any value of k. So that ([13| can be rewritten - with 
the replacement X''f{X) f{X) - like 

which shows explicitly a zero of order k + 1 a,t the 
origin, as mentioned earlier. 

Finally, under the scale transformation X aX, 
the Lagrange formula writes, after the change of vari- 
able t ^ at: 



FiaX) 



X 
i^k\ 



jyj{a-tfd]t'[X'f>{Xt)] 



(15) 



B. Evaluation of amplitudes with test functions 

From now on, we shall work with super regular test 
functions chosen among the partitions of unity. The 
property that the integral ^ is independent of the 
precise choice of this partition of unity will come out 
explicitly [T3j[14]. Moreover, if / is a PU with a given 
support, any power of /, /", is also a PU with the 
same support. They are said to be equivalent. In 
the limit where f 1 over the whole space, they 
are strictly equal. This property is essential in order 
to have in ([2]) a field independent of the construction 
of / and preserve the basic Poincare and Lorentz in- 
variances [12]. One possible realization of a SRTF 
partition of unity is detailed in Appendix [Xj By con- 
struction, the test function is 1 everywhere except in 
the vicinity of the boundaries. 

To begin with, we shall concentrate in this section 
on distributions singular in the UV domain. The up- 
per boundary - defining the extension of the support 
of / - is taken to be H in the UV domain, so that 
X < H. We shall denote by f{X; H) the test func- 
tion in order to keep track of this boundary. 

We suppose, in a first step, that the UV divergence 
is of a logarithmic type. This corresponds to the low- 
est order singularity with k = Q. Our formalism will 
be easily generalized to distributions with higher sin- 
gularities. Using (10 1, we have 



r°° r°° df 

A=-J^ dXT{X)X jdx[f>{Xt;H)]. (16) 

The integration limits are determined by the support 
of the test function. We thus have 



A- 



H 



dXT{X)X jdx[r{Xt;H)] . (17) 



C Test functions with running support 
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By integrating (17 1 by part on X, we have 



A 



-XT{X) 



^ dt 



f>{Xt;H) 



X=H 



X=0 



H 



dXdx [XT{X)] 



^ dt 



f>{Xt-H) 



X 



t 



(18) 



Because of the boundary conditions of the test func- 
tion, the first and third terms in the right-hand side 
of Eq.(18l are equal to zero, and it remains 



H 

A = I dX dx [XT{X)] 

'0 



^ dt 
~t 



r{Xt-H) . (19) 



Up to now, we have only made use of identities, and 
of the fact that the test function should be zero at 
some boundary defined by Xmax = H. Calculating 



( 19 ) with a constant H leads to the usual result using 
ordinary cut-ofi^, i.e. an amplitude which behaves like 
Log[H]. The limit / ^ 1 would thus be achieved by 
letting — > oo, leading to a logarithmic divergence. 
This is equivalent to the standard procedure using a 
sharp cut-off on the variable X. Clearly one should 
consider another construction of the test function in 
the UV - as well as IR - keeping track, in the ultimate 
limiting process, of the generic scaling properties em- 
bodied in any PU. 



C. Test functions with running support 

In order to go beyond this naive, infinite if- limit, 
we shall consider a boundary condition for which the 
boundary H depends on the running variable X |12] . 
We thus define a function g{X) by: 



H{X) = TfXg{X) 



(20) 



up to an additive arbitrary finite constant irrelevant 
in the UV domain. This function depends on an ar- 
bitrary dimensionless scale rj^ which is related to the 
shape of the test function near H . In the next subsec- 
tion we detail the conditions to be imposed on g{X). 

The extension of the procedure explained in |IIB| 
to deal with running test functions should be done 
with care. Replacing directly H by H{X) and 
f^{X\H) by a single function depending on X 
with F>{X) = f>{X; H{X)) would not make any dif- 
ference since F^{X) is a function which vanishes at 



the upper boundary denoted by X^ax and thus ( 19 ) 
will still behave like Log[Xmax], which is divergent 
when Xmax — > oo. 



Clearly, Lagrange's formula in (12 1 applies for a 
given X-value at a fixed support. So that in the t- 
integral, Xt just moves the argument of / along this 
fixed support up to a value of t such that Xt — H{X). 
In order to keep track of this fixed support, it is conve- 
nient to consider the function _F^of the two variables 
X and Y defined by 



F>{X,Y) = f>{X-H{Y)) 
and the test function / is now given by 
f>{X;H{X))= F>{X,Y) 



(21) 



(22) 



Each specific Lagrange formula is then applied to the 
X dependence of F only. In the UV domain for in- 
stance, it leads to: 



F>{X,Y)L 



Y=X 



dtdt [F>{Xt,Y)\y^^ 



-X 



dt 

T 



dx[r{Xt,H{X))] 



^XH'{X) 



jdYF>{Xt,Y)\^^^. (23) 



We shall construct the test function in such a way that 
the last term on the right hand side of ( 23 1 vanishes 



because in the UV there will be no overlapping of the 
domain where dyF^ {Xt, Y) is finite with the domain 
of integration on t. 



D. Extension of singular distributions 

1. In the ultraviolet domain 
Following the general procedure detailed in Sec. 



II B we can write, with (20 1,(21 1 and (23 1 



dX dx [XT{X)] / ^F> {Xt, X) . (24) 



In order to extend the test function to 1 on the whole 
space, we shall consider a set of function g{X), de- 
noted by ga{X), where by construction a is a real 
positive number less than 1, and where the limit 
a \^ corresponds to / ^ 1. We note that, since 
dx [XT{X)] will give by itself a finite X-integral when 
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Xmax oo (the singularity of T-^{X) has been cho- 
sen of order fc = 0), the process of extending the test 
function to unity over the whole integration domain 
in X may be chosen most simply to result, first, in an 
X-independent t-integral and, second, to X^ax — > oo, 
where X„iax — H{Xmax), and vice versa. 

Different choices of ga{X) are possible, provided the 
above two properties are satisfied. Anticipating the 
discussion of the next subsection, it should also pro- 
vide a meaningful extension of singular distributions 
in the IR domain, as detailed in Appendix |A 2[ Fol- 
lowing |12] we shall choose, as an example 



-1) 



(25) 



One might also consider other choices, like ga{X) — 
\^{a—X)X or Exp[(a — The difference between 
these choices will mainly be in the rate at which Xmax 
goes to infinity when a approaches 1~, the extended 
distribution being the same, as it should. 

To see that the limit a ^ 1~ corresponds also to the 
extension of the PU to the whole integration domain, 
it is sufficient to look at the maximal value of X at 
which the PU test function goes to zero. It is defined 
by X„^ax = H{X^ax) so that, with ([25 1 



(26) 



The limit X^ax ^ oo is achieved by the condition 
a ^ 1^, with > 1. The running support of the PU 
test function stretches then over the whole integration 
domain. 

For finite values of X, lim^^i- ggjX) = 1, so that 



the upper limit of integration on i in (24 1 is just rj^. 
Near X^axj i-G. in the asymptotic region where the 
test function is less than 1 and goes to zero ( cf Ap- 
pendix A2 and Fig. 6), we have ga{X) ~ This 
region of integration, however, does not contribute at 
all to the amplitude A since the integral on X is now 
completely finite and thus insensitive to the limit of 
very large X, when X c± Xm ax oo. 

V 



We can thus take in (24 1 the limit a 
physical amplitude writes 



The 



A 



dt 
T 



dX T>{X) 



with 



T>{X) = dx[XT{X)]Log{v^) . 



(27) 



(28) 



The extension T>{X) of the distribution T{X) in the 
UV domain depends logarithmically on t he a rbitrary 
scale T]^, with rj^ > 1. The amplitude (27 1 is now 
completely finite. 

The generalization of this procedure to singular dis- 
tributions of any order can easily be done in a very 
similar way [12^. It leads to the following extension 



T>{X) 



i-xl 

kl 



a^+i [xnx)] 



dt 
~t 



(l-^)^ (29) 



for /c > 0, while T>{X) = T{X) for fc < 0. 

Note that we do not need the explicit form of the 
test function in the derivation of the extended dis- 
tribution T^{X). We only rely on its mathematical 
properties and on its running construction. 



2. In the infrared domain 

The extension of singular distributions in the IR 
domain has been studied in details in Ref. [Sj. We 
recall here how it should be understood in terms of 
SRTFs, as proposed in [12]. In order to keep track 
of the explicit regularization of the X = Q singularity 
by the test function, we can study the amplitude in 
terms of [/^]^ ^ and write 



A^ / dXT{X)[J<{X)f 
lo 



(30) 



We can thus apply the Lagrange formula (13 1 to one 
of the /^'s only, and get 



A 



dXT{X)f<{X) 
{1-tr 



X^ 

kl 



dt- 



(31) 



By a change of variable Xt ^ X we have 



A 



r 



dX- 
X 

T 





'X' 


[ dtT 




'q 


T 



9^+1 F<[X,Y]\ 



y=4 



(32) 



Using the boundary condition on resulting from 
the behavior of the test function near the origin and 
detailed in Appendix |A2[ we can write, after integra- 
tion by part and in the limit where the test function 
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extends to 1 on the whole space "'^ 



A = hm / dX T{X) J{X) 



(33) 



with 

f<{X) = {-lf+'d'l^] 
1 



k\ 

t''+^ \ t 



(34) 



Note that the derivatives in ( 34 1 and ( 35 1 below have 



to be taken in the sense of distributions (c/ Appendix 
B3). The scale fj is positive and arbitrary. 

For an homogeneous distribution, with T[X/t] = 
t''T{X), the t integration can be carried out to give 



T<ix) = (-i)'=a^+^ 



-TiX)Log{ijX) 



(_^ \k 

^—^HkC^5^^\X) 



(35) 



with 



T{X)X''dS 



(X=l) 



7- 



and ■0 is the usual digamma function with ■0(1) 

The extension T'^{X) differs from the original dis- 
tribution T{X) only at the singularity. 



III. APPLICATION TO LIGHT-FRONT 
DYNAMICS 

A. Covariant formulation of light-front dynamics 

In CLFD, the state vector of a physical system is 
defined on the LF plane of general orientation to-x = ^, 
where lu is an arbitrary lightlike four-vector and 
^ is the LF " time" [11^ . Standard light- front dynamics 
is recovered by choosing u — (1,0,0,-1). We shall 
take C = 0; for convenience. 



Any bound system is entirely described by its state 
vector (p) . It corresponds to definite values for the 
mass M, the four-momentum p, and the total angular 
momentum J with projection a onto the z axis in the 
rest frame, and forms a representation of the Poincare 
group. It depends on the position lu of the light-front. 
The four-dimensional angular momentum operator J 
is represented as a sum of the free and interaction 
parts: 

Jpi^ — j'fj + Jpu ■ (36) 
In terms of the interaction Hamiltonian, we have 



jint 

pu ~ 



W\x){xpUJy - Xyijjp)5{uj-x) d^x . (37) 



From the general transformation properties of both 
the state vector and the LF plane, it follows [TOj fTT] 
that 



where 



(38) 



(39) 



The Eq. (38) is called the angular condition. 



This equation does not contain the interaction 
Hamiltonian, once (p satisfies the Poincare group equa- 
tions. The construction of the wave functions of states 
with definite total angular momentum becomes there- 
fore a purely kinematical problem. The dynamical de- 
pendence of the wave functions on the LF plane ori- 
entation now turns into their explicit dependence on 
the four-vector lo. Such a separation, in a covariant 
way, of kinematical and dynamical transformations is 
a definite advantage of CLFD as compared to stan- 
dard LFD on the plane t + z = 0. 

According to the general properties of LFD, we de- 
compose the state vector of a physical system in Fock 
sectors. Schematically, we have 



ID 



|2) 



(40) 



Each term on the right-hand side denotes a state with 
a fixed number of particles. In the Yukawa model, the 
analytical form of the Fock decomposition is 
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"(P) = Er^^E / 0n,.a'(fci...A:„,p,c.r„)a^,(ki)c^(k2)...ct(k„)|O) 

n=l a' 



5^'^\ki + . . . + fc„ - p - UITn)2{uj-p)dTn ]^ 



(41) 



where 4'n,aa'{---) is the rt-body LF wave function, 
the so-called Fock component, describing the state 
made of one free fermion and {n — 1) free bosons, 
(c^) are the free fermion (boson) creation operators, 
Efe, = y^kf+^nf , and m/ is the mass of the parti- 
cle I with the four-momentum fc;. The combinatorial 
factor l/(n — 1)! is introduced in order to take into 
account the identity of bosons. We do not consider in 
this approximation excitation of fermion-antifermion 
states ("quenched approximation"). 

The variable r„ describes how far off the energy 
shell the constituents are. In practical calculations, 
the infinite sum over n is truncated by retaining terms 
with n which does not exceed a given number TV, 
while those with n > N are neglected. Decomposi- 
tions analogous to Eq. (41 1 can be easily written for 



QED [18j or for a purely scalar system [T5] . 

The construction of the spin structure of the wave 
functions 4)n,aa' is very simple, since it is purely kine- 
matical. This structure should incorporate however 
w-dependent components in order to fulfill the an- 
gular condition ( 38 1 . It is convenient to decompose 
each wave function 4>n,aa' into invariant amplitudes 
constructed from the particle four-momenta (includ- 
ing the four- vector w!) and spin structures (matrices, 
bispinors, etc.). In the Yukawa model we have for 
instance, for N = 2, 



= (pi u„-{ki)uc,{p) 



(P2 + f2 



uj-p 



(42a) 
uAp) , (42b) 



since no other independent spin structures can be con- 
structed. Here it's are free bispinors of constituent 
mass m; ipi, tp2, and are scalar functions deter- 
mined by the dynamics. 

The eigenvalue equations for the Fock components 
can be obtained from the Poincare group equation 



(43) 



by substituting there the Fock decomposition ( 41 1 of 



the state vector (j){p) (here and below we will omit, 



for shortness, all indices in the notation of the state 
vector) and calculating the matrix elements of the op- 
erator p2 in Fock space. After decomposition of the 
momentum operator in free and interacting part, we 
can easily get the eigenstate equation |T5]: 



(44) 

where iJ"** is the interaction Hamiltonian in momen- 
tum space: 

if™*(wr) = / iJ"*(x)e-'("'^)^rf''a; . (45) 



According to the decomposition (41 1, the conservation 



law for the momenta in each Fock component has the 
form 



kn= P+ l^Tn 



(46) 



Hence, the action of the operator 

i^pmj -M^onthe 

state vector reduces to the multiplication of each Fock 
component by the factor — M'^ — 2{uj-p)Tn. 

It is therefore convenient to introduce the notation 



g{p) = 2{uj-p)f(t>{p) 



(47) 



where f is the operator which, acting on a given com- 
ponent 4)n,aa' of 4>{p) , givCS Tn4>n,acr' ■ G (p) haS the 

Fock decomposition which is obtained from Eq. (41) 
by the replacement of the wave functions 4'n,cTa' by 
the vertex functions r„ (which we will also refer to as 
the Fock components) defined by 

Ua'{ki)TnU„{p) = (s„ - M^)(/)„^o.^/ , (48) 

with s„ — {ki + . . . fc„)2. The vertex function r„ is 
represented graphically by the n-body function shown 
on Fig.([l]). 

Since for each Fock component s„ — Af^ = 2{uj-p)Tn, 
we can cast the eigenstate equation in the form 



Gip) 



1 

2^ 



dr 

-H"'\u;t) —Gip) . (49) 
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1) bosons 



FIG. 1: Vertex function of order n. 



The physical bound state mass M is found from the 
condition that the eigenvalue is 1. This equation is 
quite general and equivalent to the eigenstate equation 
P'^cj) = M'^cj). It is nonperturbative in the sense that it 
sums up to all orders irreducible contributions which 
involve at most TV particles at any given (light-front) 
time. 



B. Eigenvalue equation 

According to our discussion in Sec. a test func- 
tion should be attached to each particle creation or 
destruction operator in momentum space. It is then 
easy to extend the eigenvalue equation we developed 
in [THlIin] in order to account for these test functions. 
We can therefore apply the following rules: 

• To each external boson or fermion line of mo- 
mentum ki, one should attach a factor /(kf/A^). 

• To each internal propagator with momentum kj , 
one should attach a factor [/(k^/A^)]^ . 

Since each vertex function r„ is attached to one 
fermion and n—l boson lines (external or internal) , it 
will be multipHed at least by /(kf/A^) . . . /(k^/A^). 
We can thus redefine r„ to include implicitly these 
test functions. We shall call r„ these new vertex func- 
tions: 

f „[fci . . . fc„] = r„[fci . . . A:„]/(k2/A2) . . . /(k2/A2) . 

(50) 

In the limit where all the test functions go to 1, one 
has f„ — > r„. Since the f's are SRTFs, this also 
implies that any r„ is also a super regular function 
with respect to all momenta. We will thus be able to 
apply the Lagrange formula to all loop calculations, 
and derive the extension of all singular distributions 
along the lines detailed in the previous section. 

Moreover, the only singular contributions are asso- 
ciated with the loop diagram of Fig. ([2| (or similar 




(n — 2) bosons 



FIG. 2: Typical loop diagram entering the eigenvalue 
equation for the vertex function r„. 



ones with contact interactions attached to the ver- 
tices). This generates both self-energy type diagrams 
to all orders (sometimes called rainbow diagrams) and 
overlapping loops to all orders. If one can give a math- 
ematical sense to this contribution, one can give a 
mathematical sense to the eigenvalue equation to any 
order. To show how this happens in practice, we shall 
consider in the following the Yukawa model in the 
two-body, N = 2, approximation. 



C. The Yukawa model in the iV = 2 
appr oximat ion 

The Lagrangian we start from is given by 

+ + Am'l'^' , (51) 

where Am is the mass correction given by 



Am = m — mo 



(52) 



The fermion and boson fields are denoted by ^ and <& 
respectively. In Eq. (52 1, m denotes the (observable) 
physical mass of the fermion while mo is its (asymp- 
totic) bare mass. The mass of the boson is fi. Since in 
our formalism every amplitude is finite from the very 
beginning we do not need to consider any counterterm 
to cancel divergences. We introduce however a finite 



mass correction Am in (51 1 in order to decompose our 



state vector on the basis of fields with physical mass 
m, and treat the mass term Am^'^' as an interac- 
tion. According to the general Fock sector dependent 
renormalization scheme detailed in |20j , the mass cor- 
rection, as well as the bare coupling constant, should 
depend on the Fock sector under consideration, while 
the basis on which the Fock sectors are decomposed 
will stay the same. 
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From the Lagrangian (51 1, we can easily deduce the two-body component. The functions 6i_2 are related 

to ip2 , 1^32 (|42| using ( 48 1 . The corresponding quan- 



efFective Hamiltonian on the light-front. It is given in 
Ref . [18] . It includes contact interactions arising from 
the elimination of nondynamical degrees of freedom. 
Contrary to the Pauli-Villars regularization scheme, 
these contact interactions are not cancelled by Pauli- 
Villars fermions. 

As explained also in |20j , we have to consider a pri- 
ori a specific insertion on the light-front in order to 
correct for possible violation of rotational invariance. 
This insertion is denoted by Z^. 

The eigenvalue equations are shown on Fig. ([3|. 
Since there is only one nonzero mass correction and 
bare coupling constant in this N — 2 calculation, we 
denote them simply by Am and go respectively, with 
no reference to the Fock sectors. The vertex functions 
are denoted by Fi and F2. 



titles denoted with a bar are defined according to ( 50 1 . 
The mass of the boson, of momentum k2, will be de- 
noted by /i. 

To solve the eigenvalue equation shown on Fig. 
we multiply both equations by u on the left and by 
u on the right, sum over spin indices, multiply by 1 
and i/) successively the second equation, and take the 
trace [H]. We thus get one equation /i = from the 
first equation defining Fi and two equations /2,3 = 
from the equation defining r2. Solving the last two 
equations with respect to bi and 62, we get for the 
calculation of the ground state properties for which 
M m 



61 = 2 ^1 50 ™ /(k?/A2) /(k^/A^) 



(54a) 



Ti Am Am 
=■ X < X — 



2uj.p 



Ti Am 5o 
=■ y-9- 



Ti Am 

=• K 



Ti 5o Am 



50 50 



2ijj.p 



FIG. 3: Eigenvalue equation for the Yukawa model in the 
N — 2 approximation. 



Following (42 1, we can decompose the vertex func- 



tions in the following way: 

u{pi)fiu{p) = {m'^-M^)^i 



u{ki)T2u{p) = u{ki) 



— — ui 

bi + 62 

LU-p 



(53a) 
u(p),(53b) 



where M is the bound state mass of momentum p, 
and TO is the mass of the fermionic constituent of mo- 
mentum pi in the one-body component and fci in the 



-(^1 go Am /(k?/A2) /(k^/A^ 



(54b) 



1 



90- 



1 



./(k?/A2)/(k^/A2) 



2'"' (27r)3- 

(fk'j^dx' b[{2 - x') + 2b'2{l - x') 

2x'{l — X') S' — TO^ 



where s' is the usual center of mass energy squared of 
the internal boson - fermion state: 



k', 



k', 



{1-x') 



(55) 



The usual longitudinal momentum fraction, with re- 
spect to the position to of the light- front, is given by 



(jJ-k2 

LU-p 



(56) 



The transverse momentum kj_ with respect to lo is 
defined from 



R2 = k2 



xp 



with i?2 = (i?^,ki,i?°) , (57) 



since by construction w-i?2 = 0. The last component 
in ( |57p corresponds to the longitudinal momentum R2 
with respect to co. The prime quantities are defined 
similarly. 

The vertex functions &i,2 do depend in principle on 
the kinematical variables (k_L, x) (and b'^ 2 the vari- 
ables {\i'j_,x')) of the two-body system, as defined in 
Appendix (B ll. In the limit where / — > 1 however. 



these vertex functions should tend to constants ac- 



cording to Eqs. (54). We say that they are almost 
constant. 
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The limit / 
leading to 



1 can easily be done in Eq.(54al 



bi = 2ipigom 



(58) 



On top of the three equations /i,2,3, we have to 
impose two normalization conditions The first 

one fixes the bare coupling constant as a function of 
the physical one, at the (nonphysical) on-energy shell 
point s — 



5i(s = m^) = gy 



(59) 



where Ni is the norm of the one-body state ^. This 
condition originates from the construction of the ver- 
tex function r2 in which the constituent fermion state 
(single soHd Hne in Fig. ([T|) does not correspond to 
the fully dressed one-body state but to a state dressed 
with one less boson, due to the presence of an ad- 
ditional boson in flight at the same light-front time. 
This is a direct consequence of the truncation of the 
Fock space. It should therefore be corrected by its 

/]Vi for a 



normalization, i.e. by a factor \/l — N2 
two-body Fock state truncation, where N2 is the norm 
of the two-body state. 

The second renormalization condition insures that, 
always at s = m^, i.e. on the energy shell, the two- 
body vertex function r2 should not depend on the 
position of the light-front plane, i.e. 



62(3 = m^) = 



(60) 



With this latter condition, the normalization condi- 
tion of the state vector writes EU] 



l = Ni + N2 



4mV?+4mV?ffo^2 , (61) 



with 



k^L + "I (2 - 2;') 



i\2 



k^^ -I- m'^x'^ + /^^(l — x') 



,f{kl/A')f{kl/A') . 



(62) 



It is calculated in Appendix |B 4[ This fixes ipi. The 
condition (59) with (61 1 defines go- 



90^ 9 



(63) 



^ Note a different definition of the bare coupling constant as 
compared to |20j , wliere tfie normalization correction is miss- 
ing. 



and the norms of the one- and two-body states are 

1 



N2 = 



I + 52J2 

1 + 



(64a) 
(64b) 



The condition ( 63 1 insures that the bare coupling con- 



stant is not renormalized since, in the N ~ 2 approx- 
imation, there is no radiative corrections to the ver- 
tices. The norms A^i and N2 in (64 1 are always less 



than 1 and positive, but scale dependent. Since 62 is 
almost constant and thus almost zero from ( 60 1 , we 



immediately get from (54b), 



Am 



(fk\ dx' 



{2-x') 



li'^^ + m'^x'^ + y?{l-x') 



b[. (65) 



With this value for Am, the first equation Ii = 
defines the w-dependent insertion Z^^ 
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1 



8(pim2 (27r)3 



d^k'j^dx' k^^ + m^x'{2 - x') 
(1 - x') k^L^ + m^x'^ + Ai^(l - x') 



b[. (66) 



We have identified in (65 1 and (66) with 1^9 1. Us- 



ing the general properties of bi as SRTFs, we can now 
calculate Am and as well as J2 in (62) following 



the derivation outlined in Sec.(|II|- This is done in Ap- 
pendix IS] This leads to 



Am = 



3mg2 

r-l 



mg- 



[ dx'{2-x%og 
Jo 



167r2 

m^x'"^ + ^'^{l — x') 







(67) 



(68) 



These expressions for Am and Z^ are analogous to 
the ones found in [TH]. They correspond to the usual 
calculation using PV regularization scheme - with the 

sum of different factors in taken in the limit 

Apy — > 00 replaced by a specific function in the scal- 
ing parameter rj - but without any constraint on rj 
except that it should be larger than 1. The absence of 
an w-dependent counterterm, from Z^ = 0, is a clear 
check that the TLRS does not violate rotational in- 
variance in LFD calculations, at least in the N — 2 
truncation. 
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In order to extend this approach to the more gen- 
eral case, it is necessary to calculate the contribution 
of Fig. ([2]). This is a simple extension of the calcula- 
tion of the self-energy done in this Section, apart from 
the presence of an external momentum K — X]r=3 ^« ' 
where the ki 's arc the momenta of the (n — 2) bosons 
in Fig. (H. 



IV. PERSPECTIVES 

We consider in this study a coherent framework 
based on the construction of quantum fields as op- 
erator valued distributions acting on test functions, 
focusing on the specific properties the test functions 
should obey in order to achieve a generic quantum 
field theory description which preserves Poincare and 
Lorentz invariance. These properties are of three dif- 
ferent types: 

In order to achieve the independence of any phys- 
ical amplitude on the choice of test functions, and to 
fulfill Poincare and Lorentz invariance, we choose test 
functions as partitions of unity. 

a) The generalization of the procedure to deal with 
singularity of any order, which is essential in nonper- 
turbative calculations, requires one to choose parti- 
tions of unity which are superregular, i.e. which go to 
zero with all their derivatives at the boundaries. 

Hi) Finally, the extension of distributions over 
the whole space should respect the inherent scaling 
behavior in the UV as well as IR domains. This 
is achieved by using test functions with a running 
boundary, in contrast to a fixed boundary similar 
to a naive cut-off. This scale invariance enables 
a renormalization group analysis of any physical 
observable. 

This scheme has been applied to the calculation of 
a bound fermion-boson system in the Yukawa model 
within light-front dynamics. We show, in this simple 
example, why this Taylor-Lagrange renormalization 
scheme is very natural to use in nonperturbative calcu- 
lations on the light-front. It is however quite general, 
and can be applied to any formulation of quantum 
field theory. 

The mathematical properties of the test functions 
translate immediately to similar properties for the 
many-body vertex function describing the Fock com- 



ponents of the state vector. The calculation of the 
anomalous magnetic moment of a fermion, in the 
simplest two-body truncation, can be done very eas- 
ily in the Yukawa model, following the derivation of 
Ref. [20], but without any reference to Pauli-Villars 
fields and infinite mass limit. 

Our formulation can form the basis of more involved 
calculations including higher Fock state components, 
like for instance the calculation of the anomalous mag- 
netic moment of a fermion in a nontrivial three-body 
Fock state truncation. It can also be extended to the 
investigation of baryon properties at low energies in 
chiral effective field theory on the light-front [21]. 
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APPENDIX A: CONSTRUCTION OF THE 
TEST FUNCTION 

1. Partition of unity 

A partition of unity on the interval [a, b] is built up 
from a family of functions Pi{X) with 



N 



PiiX) = 1 for any X G [a, b] . 



(Al) 



A simple realization of this condition is to choose a 
basic function u{X) such that 

u{X) + u{h-X) = l for any X € [0, h] , (A2) 

where h is an arbitrary positive real number, and take 

(3i{X) = u{\X ~ ih\) for \X-ih\<h. (A3) 

One typical example of such construction is shown 
on Figs.(4[5l. It corresponds to the following basic 
function 
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\X-ih\ 



dv exp 



for 
for 



iX-i/il < h 
\X -ih\>h , 



(A4) 



with the normahzation factor J\f given by 



J\f- 



dv exp 



(A5) 



The parameter v is an arbitrary positive real num- 



has thus the following properties 
-X) 



fpu{X) 




for X e[Q,h] 
for X e [h, Nh] 
Nh) for X e [Nh,{N + l)h] . 

(A7) 




FIG. 5: Construction of a partition of unity from the basic 
functions l3i{X) given in ( A4l, for TV = 5, = 1 and u = 1. 
The partition of unity fpuJX) is the sum of all the basic 
functions. 



ber governing the shape of the function [3i{X), from 
a triangle like shape to a square like one, with all its 
derivatives equal to zero at the boundaries given by 
X={i- l)h and X l)h. The function 



N 



fpu{X)^Y.P^^^) 



(A6) 



2. Extension with a running support 

The construction of the running test function intro- 
duced in Sec jll C| can be done using the basic functions 
presented above. All the properties of the functions 
Pi are preserved when h depends on X, with h given 

by 



h{X) = 7fXg^{X) + {a-l) 



(A8) 



The additional constant (a — 1) is chosen such that, 
in the IR domain, the lower boundary of the support 
of the test function goes to zero when a goes to 1 . 

With this choice of h{X), it is sufficient to consider 
only two functions Pi, with i = 0, 1, so that the test 
function is given by 



fiX;h)^Po{X)+(3iiX) 



(A9) 



The basic functions building up unity are shown on 

Fig.del. 



2 Extension with a running support 
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FIG. 6: Construction of a partition of unity with running 
support extension, for i = (dashed line) and i = 1 (solid 
line), with a = 0.95 and /x^ — 2. The upper curve shows 
the IR domain while the lower curve shows the UV domain. 



This construction is of course not unique but is just 
mentioned here as an example of the construction of 
a running test function. We recall that the explicit 
form of the test function / is indeed not necessary in 
the calculation of the physical amplitudes, as shown 
in Sec. |TT1 

In the UV domain, i.e. for X oo, the boundary 
of the test function, denoted by H, is deduced from 
the boundary condition on i.e. from X — h = h. 
We have thus 

H{X) = 2h{X) ^ 2rj^Xga{X) + 2{a - 1) . (AlO) 

For large X, and after redefinition of the arbitrary 



scale 2ri ^ rj , we have 



H{X) ~ r^^X" 



(All) 



which is precisely of the form given in ( 20 ) . The maxi 
mal value of X, defined by X^a 

by 



H{X,nax), is given 
(A12) 



It goes to infinity when a — > 1 , with rj^ > \. 
hmit / — > 1 is thus achieved by the limit a ^ 1" 



The 



In the IR domain for any X close to 0, the test 
function is for X < Xmino with Xmmo given by 



-^minO — 2/l(A^mirio) — 2ry -^mitiO 



+ 2(a- 1) , (A13) 



so that, for a — > 1 



X, 



minO 



2(1 - a) 
2rf 



1 



(A14) 
1 with 



In turn the test function is 1 for X > Xmm 
given by 

Xmtnl = h{X,nrnl) = Kyiinl + [a ~ I) , (Al5) 

SO that, for a — > 1^, 

(1-a) 



X 



mini 



— 1 



(A16) 



It also goes to zero like (1 — a). Hence when a ^ l~ 
the test function effectively tends to a step function. 
In this region where the test function is less than 1, 
we have (1 — q)X > /i((l — a)X). This defines the 
lower limit for the t integration in Eq. (34 1, with (1 — 
a)X/t > h{{l - a)X/t) i.e. 



t > X {if - 1) 



(A17) 



Since 77^ is an arbitrary real number larger than 1 , we 
can define = 77^ — 1, with fj an arbitrary positive real 
number. Wc have in that case t > rjX. 



APPENDIX B: CALCULATION OF THE 
VERTEX FUNCTIONS 

1. Kinematics 

It is more appropriate to calculate all momenta in 
the reference frame where p = since this reference 
frame is the same whatever the number of constituents 
we have in the Fock sector. For the calculation of the 
self-energy type diagram, we denote by ki (k2) the 
momentum of the intermediate fermion (boson). They 
satisfy the on-mass shell conditions, and their vector 
parts are decomposed in parallel fcy and perpendicular 
fcj^ components, relative to the position of the light- 
front u) = (cjojw). We have therefore, with = 



LO-p = UJ{)M , 

uj-k2 = X u>oM 



- wo(£fei - fci.ii) 



(Bla) 
(Bib) 
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with = ^/k-]^ + ^2 II + M^- We thus get 



and 



SO that 



2xM 



2 ^ ^2 ^^2 



2xM 



k\ + ^2 xM 



2xM 



(B2) 



In the hmit of large transverse momentum, or for x 
0, we have 



2 _ (ki + M 



2\2 



Ax^AP 



(B3) 



The momentum ki of the intermediate fermion can 
be calculated easily with the replacement x ^ 1 — a;, 
m fi and kj^ —^ — kj^, so that 

, , / k2 +m2 (1 - x)M\'^ 2 

k? = / , — + ^ ' — - . (B4) 

In the limit of large transverse momentum, or for x 
1, we have 



2^ (ki+m2)2 



4(1 - x)2Af2 



(B5) 



The mass correction 



The equations for the two-body components lead 



The test functions will alternatively correct the be- 
havior of the integrand of / for x' ~ 0, with f(k'2) or 
for x' ~ 1, with f(k'i). We can thus divide / in two 
contributions /1.2 for x' < 1/2 and x' > 1/2. Keep- 
ing only the test function which is not identically 1 in 
each interval, we have, for Ii for instance: 



I.^Tdx' Hd^k', /(kl/4.-^-A-) 

Jo Jo k^2 + m^x'^ + - x') 

(B9) 

where the argument of the test function is valid only 
in the h'^ ^ 00 region. With the following change of 
variable 



X = \s! \l2mAx' , 



(BIO) 



we have 
/i = 27rmA 

J 

dX 



fix') 



2mAx'X + m2x'2 -|- ^'^{l - x') 



(Bll) 



To get rid of any momentum scale in the behavior of 
the integrand, we make the change of variable 



Z = 

so that Ii writes 
A 



ni^x''^ + /i2(l — x') ' 



(B12) 



h = 27r- 



dx X 



dZ^-^ — , — . B13 



to the expression (|65| for the mass correction. Using Using ^ with (|15|) and ([U), for fc = 0, we have 



( 54a I , it is enough to consider the following integral 



/ dx' ( d'k\ -4^ 
Jo Jo k'i 



-|- m?x'^ + /i2(l — x') 

(B6) 



Ji = 2tt 



A 



^00 

x I dZdz 
'0 



2A 



'^Z+l 



where k'^ and k'2 are given in Eqs. (B2 B4), in terms 



of x' and k'^, and for the most dangerous region in 
k'^ (k'^ ^ 00) and x' {x' 0, 1) by (|b3| and ( |b5| . 

Since / has no singularities in = 0, 1, we can even 
consider: 



with 



r '^F{ZH\Z') , (B14) 
rn^x''^ H- /i2(l - a;') 



(B15) 



k1 



k' 



k'l 



4(1 -x'Y-m? ' 



4x'2t 



(B7) 
(B8) 



Because of the structure of the test function, the inte- 
gration domain in t is limited according to the support 
of F, i.e. 



ZH"^ < H{Z') « ■q'Z'gaiZ') • 



(B16) 



2 The mass correction 
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Moreover, the integrand for the Z-integration behaves This gives the result indicated in Eq. (67 1, 
now as when Z oo, i.e. the integration over 

Z is convergent. We may thus put / — > 1 (by letting 
a go to 1~) and the only remaining trace of the test 
function is in the scale 77. Doing the integration over 

Z, it remains 3_ The ci;-dependent insertion 



TT I dx' 
'0 



" dt 

T 



TT I dx' Log 

10 



■q m 



m'^x'^ + fi^{l — x') 



We can do the same calculation for the second part I2 
of I and get a similar result, the only difference being 
the integration limit on x' from 1/2 to 1. We thus get 
for /: 



TT Log[r/] — TT 



/ dx' Log 
^0 



m^x'^ + /i2(l - x') 



The equation for the one-body component leads to 
the expression (66 1 for the w-dependent insertion. It 



is enough to consider the integral 



J 



/ dx / d(k'l)/(kiVA^)/(k^VA^) 
Jo Jo 



k'_J + m'^x{2-x) 



(1 - a;)[k^^ + m2a;2 + ^^{1 - x)] 



(B17) J can be expressed as 



(B18) 



J = 



f ^ rrf(k'i)/(k',VA2)/(k',VA2)+ fdx ^d{k'l){2m'x-^,')-l 
Jo Jo Jo Jo k^ 



J1+J2 



/(k'//A2)/(k^VA2) 
^ + m'^x^ + IJ?{1 — x) 

(B19) 



Ji is well defined since the test function takes care 
explicitly of the behavior at high transverse momenta 
and at x = 1 or at a; = after the change of variable 
a; — > 1 — a;. W ith th e asy mptotic expressions at large 
k'^^ given by (B3), (B5l and after an evident change 
of variable 



Ji = 2r7iA 







f 








x"^ 





fj. -m \2 



' 2mA 

(l-a;)2 



(B20) 

The first test function treats the singularity at a; = 0. 



It is sufficient to apply the extension formula ( 34 1 for 
k=0. With T{x) = 1/x, we have 



f<ix) 



-dx ^ 
da;Log[fix] 



dt t 
¥ X 



where fj is an arbitrary (dimensionless) scale. The 
derivative of the Log function should be understood 



in the sense of distributions: 



d f°° 
< —Log{x),ip Lim^^o J Log(a;) [-ip'{x)] dx 

ip{x) 



= Lim, 



(^(e)Log(e) + 



-dx 



(B21) 



With Lp{e) = ip{0) + €ip'{^), and since ip'{^) is finite, 
we have 



Lim^^o 



fix) 



dx + v3(0)Log(e) 



(B22) 



This is precisely the definition of the pseudofunction, 
denoted by Pf , of 1/x introduced in [T3] , so that 



f<{x) = Pf 



(B23) 
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Hence J7i writes, after the change of variable ^ 



1-X ' 



Jl 



2mA I dxPfi- 

10 



— /( — ) 



,1, 



= 2toA / dx{l - x)Pf{-) {' 
Jo X 

= VLog(x)r=^ = 0, 



W 



w 



since the remaining test function defines Pf{-^) and 
its integral outside the singularity is just — The 
notation indicates simply that f{u) should 

be taken at the value u = a, the lower limit of in- 
tegration being taken care of by the definition of the 
pseudofunction . 

We can now proceed to the calculation of J2 ■ With 
the same change of variable in the k'^^-integration 
leading to ( B20 ) , we have 



J2 - 



pi poo 






dx dXf 




/ 




x^ 




2m2(l - x) - 







2mA 



(i-xy 



X 



2mA ^ 2A^(2 ^) 



dx [ — — 

Jo 

4y2A2x2 



4r2A2(l-a;)2_ 

2to2(i - - 
l + Yx{^J? ~m^{2-x))' 



The arguments of the test functions are here taken for 
Y close to since this is the domain where they differ 



from 1. 

The product of the two test functions being invari- 
ant under the change a; — > 1 — x, it can be represented 
by a single function of argument x{\—x) with the same 



support, say G 



the Lagrange formula ( 13 1 applied on G gives 



With Z{x) 



/oo 
dtdtG[Z^ix)t^] 



Integrating by part on x, J^2 finally writes 



J2 = 



dx 



'dY 
Y 

x{l ^ 



dx 



"dt 



G 



Y^x^il^xY 



v){2m^{l-x) -^2) 



_(1 - 2x){l + Yx{^? - to2(2 - x)))_ 
The test function G treats the singularity at 1" = 
by the extension of ^ to Pf{Y). Thereupon its sup- 
port can be extended to infinity, the upper integration 
limit in t is restricted to a constant rj and G becomes 
unity on every domain of integration. The remaining 
integration on x is then trivially zero. 



Normalization of the state vector 



The calculation of (62) is very similar to the cal- 
culation of the mass correction Am detailed in the 
previous subsection. We can simply rewrite J2 as 



J2 = (4m2 - n^) dx' 



d^k\ 



+ m'^x'^ + ^^2(1 — x') 



dx' 



d^k', 



/(k',VA2)/(k^VA2) 

'^^^ + m^x'^ + ^1^(1- x') 



The first term is convergent (hence the limit / ^ 1 
taken) , while the second one is exactly / calculated in 



( B 2 1 . It depends logarithmically on i]^ . Its explicit 
expression is not needed in our calculation. 
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